A strange metal with a small Fermi surface and strong collective excitations 
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We develop a theory of a hybrid state, where quasi-particles coexist with strong collective modes, 
taking as a starting point a model of infinitely many one-dimensional Mott insulators weakly coupled 
by interchain tunneling. This state exists at an intermediate temperature range and undergoes an 
antiferromagnetic phase transition at temperatures much smaller than the Mott-Hubbard gap. The 
most peculiar feature of the hybrid state is that its Fermi surface volume is unrelated to the electron 
density. We present a self-consistent derivation of the low energy effective action for our model. 
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I. INTRODUCTION 



The presence of strong collective modes interacting 
with quasi-particles is a distinctive feature of many 
strongly interacting systems such as 'bad' metals, weakly 
doped Mott insulators (such as the cuprates) and heavy 
fermion materials. This interaction is believed to result 
in a variety of unusual phenomena observed in these sys- 
tems such as the violation of the Mott-Regel limit in 
the temperature dependence of the electrical resistivity 
of bad metals or the absence of quasi-particle peaks in 
the spectral function of the cuprates. The lack of non- 
perturbative techniques in dimensions higher than one 
makes a detailed theoretical description of these phenom- 
ena quite challenging. One sucessful approach has been 
developed by Chubukov, Schmalian and Abanov, who 
have studied the so-called spin-fermion model put for- 
ward by D. Pinesi. This model is semi-phenomenological 
and postulates the existence of a strong, coherent, collec- 
tive mode, which interacts with quasi-particles located in 
the vicinity of a Fermi surface. This model is reviewed 
comprehensively in Refs 2] and has proven quite suc- 
cessful in explaining various properties of the cuprates. 
However, a derivation from a microscopic Hamiltonian is 
lacking. 

In this paper we provide a microscopic derivation of 
a model in the same class as the spin-fermion model of 
Pines and Chubukov. Namely, we continue to develop a 
theory of a hybrid state combining features of a Landau 
Fermi liquid and a Mott insulator suggested in Ref.|3]. 
This state is characterized by the coexistence of well- 
developed collective modes with quasi-particles. The lat- 
ter ones have a small Fermi surface (SFS), the volume of 
which is unrelated to the total number of electrons. By 
definition, the Fermi surface (FS) is small if its volume 
is less than the maximum volume allowed by Luttinger's 
theoremiiSiSi The existence of such a state does not 
contradict Luttinger's theorem since the latter, contrary 
to popular belief, does not fix the volume of FS. Instead 
the theorem states that the electron density n is propor- 
tional to the volume of phase space enclosed by the sur- 
face where the single electron Green's function changes 



its sign 
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When the Green's function has zeroes, the Fermi sur- 
face constitutes only a part of this surface, namely the 
one where G(0, k) — > oo. Hence Luttinger's theorem JJJ 
does not even require the existence of a Fermi surface: 
the Green's function may only have zeroes and no poles, 
as it is the case for superconductors^ and certain one- 
dimensional systems, in which the spectral gap is gener- 
ated dynamically (for the latter case a general proof is 
outlined in Ref. [8|]). 

A metallic state with a small FS would necessarily be 
associated with a Green's function that has both poles 
and zeroes at ui = 0. In our previous work [3j we sug- 
gested a model for such a state based on the quasi-one- 
dimensional Hubbard model at half filling. The trans- 
verse hopping was treated in a Random Phase approx- 
imation (RPA). In order to understand the conditions 
of stability of such an exotic metal, one has to go be- 
yond RPA, which is the main subject of the present pa- 
per. Experimental indications of the existence of SFS 
states come from ARPES measurements in underdoped 
cuprates 9 and from the Hall effect measurements in heavy 
fermion materials 10 . 

Before turning to the calculations, we shall give a qual- 
itative account of the physics we are after. Our start- 
ing point is an ensemble of uncoupled, Mott insulating 
chains. The relevant energy scale is the ID Mott gap 
m. We consider finite temperatures T such that T -C m. 
The physics is purely one dimensional. 

We then turn on a small long range interchain tunnel- 
ing with characteristic energy scale t±. Clearly, at zero 
temperature the hopping between chains will be essential 
and induce a 3D ordered state. On the other hand, in 
the window 

t± <T,i±(k) < m , (2) 

we will recover the physics of ID Mott insulating chains. 
Here £j_(k) denotes the Fourier transform of the inter- 
chain tunneling. Furthermore, as T <C m we may to a 
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good approximation work with zero temperature quanti- 
ties in many instances. 

The crucial point is that while t ± remains much smaller 
than the Mott gap m, the Fourier transform tj_(k) can 
become comparable to m in some region of the Brillouin 
zone, i.e. we may have a situation where 



(i«r« t(o) 



(3) 



In this case an interesting "hybrid" state combining 
ID with 3D features develops. In particular, the low 
energy sector corresponds to a 3D metal with a small 
Fermi surface and quasi-particles interacting with well- 
developed collective modes. The existence of the regime 
is ensured by making the interchain tunneling long 
ranged. 

The dimensional crossover from a quasi one dimen- 
sional Mott insulator to an anisotropic 3D Fermi liquid 
as a function of the strength t± of the interchain hopping 
is sketched in Fig. ^ 
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II. THE MODEL 

The model we study is the Hubbard model with a 
strongly anisotropic hopping: 

H = - l J2 [4i,l,a C n+lM + h - C - 
n,l,fj 

I \ * ,np f 

2_^i lm C n,\,a C P,m,cr ■ 

l,m,n,p,fr 

For definitencss we consider the chain direction to be z, 
so that 1 = (l x ,l y ),m — (m x ,m y ) label Hubbard chains 
and n,p label the sites along a given chain. 

As we have mentioned before, the hopping integrals 
in the transverse directions are supposed to be small in 
comparison to t. In the limit t± = and at half filling 
the model has a Mott-Hubbard gap m. We work in a 
regime where the magnitude of this gap is much smaller 
than the bandwidth W ~ At. In our previous paper 3 
the transverse hopping was treated in a Random Phase 
approximation (RPA). In order to suppress corrections to 
RPA (at least in some temperature interval) we assume 
that the transverse hopping is long ranged (see below). 
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FIG. 1: Cartoon Phase Diagram for T <C m for weakly cou- 
pled ID Mott insulators, where m is the ID Mott gap. 



The purpose of the present work is to derive an effective 
theory for the low-energy degrees of freedom in the "ID 
Mott insulator/3D Fermi liquid hybrid" regime and to 
analyze its instabilities towards 3D order at sufficiently 
low temperatures. 



A. Uncoupled, Mott insulating chains 

Let us briefly discuss the low-energy physics for uncou- 
pled chains. In order to ease notations we suppress the 
chain index (1). Keeping only low-energy modes around 
the two Fermi points ±A;_f, the electron annihilation op- 
erators are decomposed as 

ao [ex.p(ikpx) R a (x) + exjp(-ikFx) L a (x)\ 



(5) 

where ao is the lattice spacing, x — ja,Q and kp = ir /2aQ. 
The fermionic creation operators for left and right mov- 
ing Fermions are bosonized, using the following conven- 
tions 

if a 
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Rl(x) = ^e l/ ^/ 4 exp 
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where r\ a are Klein factors with {r) a , i]b} = 25 a b and where 
/]• = 1, fi = — l.The chiral boson fields (j> a and cf> a fulfil 
the following commutation relations 



b a (x),<j) a (y)] = 2ni , a = c,s. 



(7) 



In terms of the chiral fields <j) a and 4> a we define canonical 
Bose fields $ a and their dual fields a by 



(8) 



The low-energy effective Hamiltonian density for a single 
chain takes the following bosonic form 

Hs = [(^ -) 2 + ( d ^s) 2 ] - g j j, 



16tt 



[{d x Q c ) 2 + (d x <$> c f] +<?! I. (9) 
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Here I a and I a (J a and J Q ) are the chiral components 
of the SU(2) pseudospin (spin) currents 



I+ = mm ^ 

2?r 

j+ _-mm 
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(10) 



The current-current interaction in the spin sector of @ 
is marginally irrelevant and we will ignore it in what fol- 
lows. We note that doing so enhances the symmetry in 
the spin sector from SU (2) (spin rotational symmetry) 
to SU(2) x SU(2) (rotational symmetry of the left and 
right sectors). 



1. Single- Particle Green's Function 

The single-particle Green's function for the half-filled 
Hubbard model was obtained in the framework of the 
formfactor approach in Refs In particular, when 

the charge and spin velocities are equal we have 



G (u,±k F + q) = 



Zn 



ilu =F vq 



m 



m 2 + uj 2 + (vq) 2 



(11) 



where Zq w 0.921862. In order to obtain the above ex- 
pression for Go we took into account only processes in- 
volving the emission of a single massive holon and a cas- 
cade of gapless spinons. 



2. Spin Sector 

The spin operators S" = \c i n a u Zb c n,b are expressed in 
terms of the left and right moving Fermi fields by 



si 



(-l) J 'f[4(z) ^ b L b (x)+h.c. 



[Rl(x) < 6 R b (x) + R^L] 



a [(-iyn a (x) + J a (x)] 



(12) 



The bosonized expressions for the staggered components 
of the spin operators are 



Rl(x) u a ab L b (x) 



— L=expQ$ c j tr(,/rr",. (18) 



where we have replaced the product of Klein factors by 
their expectation value 



(vivi) 



(14) 



and where the matrix field g is expressed in terms of the 
spin boson <& s and its dual field <d s by 



fa^f exp(|$ s ) iexp(-j6 s ) 
9 V 2 \iexp(ie s ) exp(-|$ s ) 
At T — we have 



(g a p(T,x)g! yS (0,0)) = 5 aS 5fa^-j=, 



«o 



(15) 
(16) 



Using (|15|l one can easily calculate multi-point correla- 
tion functions of g. 

The action describing the collective spin excita- 
tions on each chain is equivalent to the SUi(2) Wess- 
Zumino-Novikov-Witten (WZNW) model once we drop 
the marginally irrelevant interaction of spin currents. 
The WZWN action for the matrix field g(r, x) is given 
by 



W[g] = ^-Jd 2 x Tr [&>g%g] 



+ v£ j B * x Tr [9^99^gg^g] , (17) 

where x\ — vt, Xi = x, <9 M — -^jr, B is a three dimen- 
sional half-space (x% < 0) whose boundary at X3 = 
conincides with the two-dimensional (vt, ir)-plane and g 
is an arbitrary extrapolation of the function defined on 
the two-dimensional space X3 — 0, which approaches 1 at 
X3 — > —00. The action (|17fl is invariant under SU(2) x 
SU{2) transformations g -> UgW. The marginally irrel- 
evant interactions of spin currents breaks this symmetry 
down to the diagonal spin-rotational SU(2) g — > UgU^ . 
The form l|17|l of the action for the spin degrees of free- 
dom is significantly more complicated than the free boson 
representation ljUJi. The latter is very convenient for cal- 
culations in one dimension, but may be less useful when 
one considers interchain coupling due to the fact that the 
dual field Q s is non-local with respect to $ s . The formu- 
lation in terms of the matrix field g has the advantage 
of the fundamental field being the order parameter itself. 
In fact, W[g] is the Ginzburg-Landau action for a ID 
spin- 1/2 antiferromagnet . 



3. Three Point Function 

An important ingredient in our analysis are three-point 
functions of the form (Tr[g(z)cr Q ] R£(zi) L b (z 2 )). The 
large-distance asymptotics of these correlators can be 
evaluated by using the results ofi^ 



4 



(Tr[g(z)a+] R ( Zl ) L T (z 2 )> 



■ (vim) 



2ir 

■ (mm) 



(Tr[g(z)a+] e" 



(Tr[g(z)a+] e" 



\4>s(Zl) P i0 3 (z2)\ /„|0 c (Zl) „i0c(z 2 )\ 



(-x) e iM**)) H W^ e af Ko (mr 12 ) 
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Z^yjt) ^ (mr 12 )(Tr[g(z)a+] Tr[g(z + )<r-]) s , 



where zi. 2 = (ti, 2 ,:ei, 2 ), n 2 = |zi - z 2 | and z+ = (— j 
of the single-particle Green's function byi^ 



(18) 



11+12 £i±s2.-j_ ^ ne cons t an t Z is related to the normalisation Z 




(19) 



The calculation we have just carried out can be summa- 
rized by the following approximate relations 



i?t(zi)< b £ b (z 2 ) 
it( Zl )< b i? b (z 2 ) 



iZ K {mr 12 ) Tr[g(z+)<r a } 
iZ K (mr 12 ) Tr[.g(z + ) ( r a ] 



(20) 



The approximation (|2U[I fails at small distances. In order 
to remove the logarithmic singularity of Kq one needs to 
include terms corresponding to the multiple production 
of solitons and antisolitons. At energies much smaller 
than the Mott gap, the fusion l|2U|) gives rise to the spin- 
fermion vertex depicted on Fig. [21 



FIG. 2: The fermion-spinon interaction generated by fusion 



B. Long range interchain hopping 

In order to have a small parameter in our theory we 
consider the interchain hopping to be long-ranged, such 
that the Fourier transform of the hopping matrix ele- 
ments strongly depends on the wave vector. This is a 
well-known trick (see, for example Ref. an d results 

in a controlled "loop" expansion, where every integra- 
tion over the transverse momenta leads to a small factor 
Kq in three dimensions, where kq is the inverse range of 
the interchain tunnelling. The interchain hopping may 
be taken long ranged both along and perpendicular to 
the chain direction. In order to simplify the calcula- 
tions, we will constrain our discussion to the case where 
^im = tQ- ~ m )^n,p: i- e - the interchain hopping has no 
component along the chain direction and depends only 
on the distance between chains. The Fourier transform 



of the interchain tunneling is then given by 

i±(k±) = ^2h m exp(zk_L ■ (1 - m)a±). (21) 

m 

In the following we choose the interchain hopping such 
that it respects the particle-hole symmetry 



t-n,l,er < * \ 1 ) L rU,CT > 

which implies that 

?x(k ± + Q) = -*j_(ki), 

where 
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Q= — ,— 

a± aj_ 



(22) 



(23) 



(24) 



is the antiferromagnetic wave vector in the direction 
transverse to the chains. It is straightforward to general- 
ize our following analysis to non particle-hole symmetric 
cases. The basic assumptions underlying our model are 
then summarized in the following inequalities: 



W > m ~ MO) I = |f x (Q)| » *l(pj 



(25) 



Here W — At and m are the band width and Mott 
gap for uncoupled chains respectively and |p_i_a_i_|, |(p_i_ — 
Q)a_i_| S> kq. The small parameter kq characterizes the 
support of t±(k±) in momentum space. The precise form 
of the momentum dependence of t± is supposedly unim- 
portant, but in order to simplify the concrete calculations 
we shall use the following model: 



ij.(ki) 



to 



l + \h±a ± \ 2 Ko 2 ' 



\k±a±\ < 1. (26) 



Within the model i|2ti[l the integration over the transverse 
wave vectors may be replaced by integration over t = 
t±(k ± ) 



Air 2 



-fit) 



47T 



to dt 



[/(t) + /(-*)]. (27) 



Some readers may find that our approach is similar to Dy- 
namical Mean Field theory in an infinitely dimensional 
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space. This is not the case: the difference comes from the 
fact that in our model the transverse density of states is 
constant on the zone boundary. This feature strengthens 
the influence of fermionic coherent modes and utterly 
changes the physics (see the discussion in the Conclu- 
sions). 



III. PERTURBATION THEORY IN THE 
INTERCHAIN TUNNELING 



As we have already mentioned, the perturbation theory 
in the interchain tunneling can be reorganized in terms of 
a loop expansion. Every integration over the transverse 
momenta generates a small factor Kq. We will refer to 
the leading order 0(kq) in this expansion as "Random 
Phase Approximation" (RPA). 

The RPA expression for the single particle Green's 
function G was derived in Ref. Q and is given by 



G{u,±k F + q,k±) = 



G (u),±k F + q) 



1-ti 



(28) 



_(k±)G (u,±k F + q) 

Here Go is the single-particle Green's function for an in- 
dividual chain . In a purely one-dimensional Mott in- 
sulator the electron is a composite particle and as a result 
the spectral function is incoherent. Coherent electronic 
excitations reappear as soon as the interchain tunneling is 
turned on. They can be understood as antiholon-spinon 
bound states and occur at energies below the Mott gap. 
When t± exceeds a certain critical value, the dispersion 
of the coherent mode crosses the chemical potential and 
a Fermi surface is formed. As a consequence of particle- 
hole symmetry, at half-filling this Fermi surface consists 
of electron- and hole-like pockets of equal volume. A 
sketch of such a Fermi surface is shown in Fig. |3| A con- 
venient measure for the strength of the interchain cou- 
pling is given by the quantity 

Z ee M , (29) 
m 

where t is defined in - The RPA form f° r the 
Green's function features a pole corresponding to a coher- 
ent quasi-particle mode. This mode crosses the chemical 
potential when Z exceeds the critical value 

Z c = 3.33019 . . . , (30) 

and a Fermi surface is present for all Z > Z c . 

Having in hand the expression for the chain single- 
particle Green's function, we may use it to define a 
dressed interchain hopping Tr£(cj, q, kj_) by summing 
the diagrams shown in Fig. ^ This results in 



t]_(kj_) 



(31) 



l-t(kj_) G a (iJ,±k F + q) 

where the -I- sign corresponds to R and the — sign to L 
respectively. We note that the dressed interchain hopping 
is equal to the propagator of the Hubbard-Stratonovich 
field that can be introduced to decouple the interchain 
hopping 11 . 
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FIG. 3: The Brillouin zone with the electron (red ovals) and 
hole-like (blue semi-ovals) Fermi pockets of a two-dimensional 
lattice. The noninteracting Fermi surface is shown as a dashed 
line. 
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FIG. 4: The dressed interchain hopping. 



A. The Spin Sector in RPA 

In the RPA, the spin sector remains one-dimensional 
and critical. This can be seen as follows. Let us consider 
the real-space correlator between the staggered compo- 
nents of spins on different chains 1 and m 



KtiW ^ m (o)). 



(32) 



Within perturbation theory in the interchain hopping, we 
need at least one right moving and left moving fermion 
operator each on chains / and m in order to obtain 
a nonzero expectation value in the spin sector. The 
only ways to achieve this are shown in FigO Here 
the dashed lines denote the bare interchain hopping, el- 
lipses enclosing two black (white) circles represent the 
purely ID Green's function of right (left) moving elec- 
trons on a given chain and the ellipses enclosing two cir- 
cles and a hexagon stand for the three-point function 
l)18|l. Clearly all such diagrams involve at least one in- 
tegration over the transverse momentum. Hence, within 
the RPA spin-spin correlation functions remain entirely 
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one-dimensional and spins on different chains remain un- 
corrected. 




FIG. 5: Real-space diagrams that contribute to the two-point 
function of staggered magnetizations between chains (I) and 
(m). 



verge at small temperatures. Therefore RPA works only 
at finite temperatures and for its consistency the tran- 
sition temperature (below which the system is three- 
dimensionally ordered) must be much smaller than the 
Mott-Hubbard gap m. ft is therefore important to es- 
timate the corresponding corrections to RPA and their 
temperature dependence. As the first step in taking into 
account corrections to RPA we have to estimate the in- 
terchain RKKY interaction. As we have shown in section 
3-point, there is a three-point "vertex" that couples the 
spin degrees of freedom to the fermionic quasi-particles. 
In second order perturbation theory in this interaction, 
an interchain exchange interaction between the spin de- 
grees of freedom is generated. The corresponding action 
is given by 



B. Excitation Spectrum in RPA 

From the discussion above it is clear that for suffi- 
ciently strong interchain hopping Z > Z c the RPA leads 
to two types of gapless excitations 

• Fermionic particle and hole excitations over the 
Fermi surface with anisotropic 3D dispersions. 

• Collective excitations of the spin degrees of free- 
dom. These are of a purely ID nature and do not 
have a dispersion in the direction transverse to the 
chains. 

If one goes beyond the RPA, interactions between these 
two types of excitations will be generated. In the fol- 
lowing we determine the form of these interactions and 
study their effects. To go beyond the RPA in princi- 
ple requires the knowledge of the two-particle Green's 
function of uncoupled Mott-insulating chains. However, 
if one restricts one's attention to the regime of energies 
small compared to the Mott gap, the three-point func- 
tion 3point (which corresponds to a particular limit of 
the two-particle Green's function) suffices. 



S xc = Y\ dT i dx i ^2 Jlm ( Xl ~ X2 ' Tl ~ T2 ) 

x Tr[ gl (T 1)Xl ) gL(r 2 ,x 2 )] . (33) 

The Fourier transform of the leading order (in kq) ex- 
change matrix element is given by the "bubble" diagram 
shown in Fig. El where the doubles lines are the dressed 
interchain hoppings for left and right moving fermions 
and the squares denote the elements of the matrix field. 
The result is 



q+q/2 




•A '/ 
\\ '■ 
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q'-q/2 



IV. INTERCHAIN EXCHANGE AND 
ESTIMATE OF THE TRANSITION 

TEMPERATURE FIG. 6: Leading order 0(k%) contribution to the interchain 



Although it is obvious that corrections to RPA are 
of higher order in the small parameter kq, they will di- 



exchange. 



J(v,q, qj.) 



z 2 



cLu'dq' 



m 
2' 



{vq'Y 



T, 



r>9 



2' 



(34) 
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To make the calculations easier, we use the fc-dependence 
(|26|l so that the sum over the transverse wave vectors is 
replaced by integration over t according to 1)27(1 . Since 
t(k±) is peaked near zero and Q, there are two interesting 
wave vectors: qj_ = and = Q. 



A. Case 1: Z < Z c 

In this case the coherent electron modes still have a gap 
and no Fermi surface is formed in the RPA. Using l)27|) 
to carry out the summation over transverse momenta we 
obtain 



x „ f°° , arctan(ZG(s)) 
where Z = Zoto/m is defined in i|29|) and 



and 



G{s) 



Z 



ijiZq 7r 3 ao 



= a Ci, (35) 



(36) 



(37) 



As expected, this interaction is of order of k§£o- The 
numerical factor cto ranges between for Z = and 2.81 
for Z — > 3.33019. The exchange at momentum transfer 
Q is 



.7(0,0, Q) 



/ ds s 



In 



l+ZG(s 



l-ZG(s) 



(1 



=2^2 



G(s) 



-olqCi 



(38) 



where ctQ varies between for for Z = and 3.07 for 
Z^Z r .. 



Let us now turn to the exchange at wavenumber Q. We 
find that there is a logarithmic divergence in (|34|l , which 
is related to the Fermi surface formation and will be dis- 
cussed later in more detail. We regularize the divergence 
by temperature. This may be done by replacing the u/- 
integral in (|34l) by a sum over Matsubara frequencies and 
substituting the finite-temperature Green's function for 
their T = analogs in the dressed interchain hoppings 
Tr.l in l|34|) . At low temperatures the single particle 
Green's function is given by- 



G R (uj n ,q) = 



-oo 



47T 2 7 



xRe 



2iB - - - 



\U3 - vq\ 



1 2 2 

i w — v q 



4 2 2nT\uo -vq\ '2 



(41) 



The singular piece of J slng (0,0, Q) diverges logarithmi- 
cally with temperature and is estimated as 



J sing (0,0,Q) m -Ci In 



s+ ds 



(i 



where s± are solutions to the equation 
\-ZG(s±) = . 



(42) 



(43) 



In order to establish the exchange at nonzero values of u 
and q we have calculated J(u>, q, Q) numerically at small 
temperatures. Rather than using the finite-T Green's 
function l|41|l we work with the T — expression i|ll|) 
and replace u by the discrete Matsubara frequencies. For 
small temperatures this is a reasonable approximation. 
We find that | J(u, q, Q)| is largest at u = = q. 

In addition to the singular piece l(42|l there also is a 
regular contribution to the exchange. As long as we are 
close to the transition, i.e. 



B. Case 2: Z > Z c 

In this case a Fermi surface in the form of electron 
and hole pockets is present. The presence of zero en- 
ergy modes does not really affect the exchange at zero 
momentum transfer, which is given by 



.7(0,0,0) wCi / ds 



sf(s) 



G(s) 



a' d , (39) 



where 



f(s) = 2arctan(£(s)ZG(s)) - arctan(ZG(s)) , 
£( S ) = min(l, [ZG(s)}- 1 ) . (40) 

We find that a' Q starts at 1.405 for Z -> 3.33019, then 
goes through a maximum of approximately 1.48 around 
Z w 4.18 and then diminishes slowly. Hence the ex- 
change at wave number zero play a subdominant role. 



< 1 



(44) 



we may estimate the regular contribution to the exchange 
by its value at the critical strength Z c of the interchain 
tunneling. The latter is given by (|38H 



J rcg (0,0,Q) sa -3.071 d. 
In the regime (|44|l 

J sing (0, 0, Q) w -1.318 ^Z-Z c Ci In 



(45) 



(46) 



The total exchange constant at wave number Q is then 
estimated as 



J(0, 0, Q) w -C x 3.071 + 1.318 y/Z - Z c In 



T 



(47) 



Having determined the exchange constant we are now 
in a position to estimate the temperature at which a mag- 
netic instability develops. In the absence of interchain 
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hopping the correlation functions of the matrix field at 
T = are given by JTJJ. At T > we have 



(9a0(T,x)g+ s (O,O)) = SagSpy 



21 sinh i 



■ IVT 



]) 



r,(48) 



if we neglect the marginally irrelevant current-current in- 
teraction in the spin-sector of the Hamiltonian (JSJ de- 
scribing the ID Mott insulating chains. If one takes it 
into account in renormalization group improved pertur- 
bation theory one obtains^ 



(tv(ga a ) tr(ffV)) = 5, 



a/3: 



In [A] ttTclqv- 1 



sinh (^[a; + ivr]) | ' 



(49) 



where A is a high-energy cutoff, which we may take to 
be of the order of the hopping integral along the chain 
direction. Carrying out an analogous calculation for the 
dimerization operator we find 



(tr(p) tr(fft)) = 5 afl 



[in [A]] " TrTa v- 



sinh(^ 



(50) 



Upon Fourier transformation and analytical contibuation 
one finds the following result for the dynamical magnetic 
susceptibility of the uncoupled chain system^ 



Xid(oJ,q) 



a 



:uj±vq] 
1 AttT J 



2T 



r[ 



4ttT 



Til 



3 _ 

4 ' 4ttT J 



(51) 



The dimerization susceptibility is equal to l|51|) apart 
from the prefactor, in which v/rn[A/T] is replaced by 

[In [A/T]] _3/ ' 2 . Hence the staggered spin susceptibility is 
always more singular than the dimerization susceptibil- 
ity and as a result the dominant instability of the spin 
sector is towards Neel order. The enhancement of the 
spin susceptibility as compared to the dimerization sus- 
ceptibility is caused by the leading irrelevant operator in 
the Hamiltonian, namely the interaction of spin currents. 
If we were to add an interaction to the underlying lat- 
tice Hamiltonian in order to eliminate this interaction, 
the symmetry between dimerization and the staggered 
components of the spins would be broken by some other 
irrelevant operator. The the dynamical susceptibility of 
the coupled chains system can be determined by a ex- 
pansion in the interchain coupling of the type discussed 
iniSii 3 .. The leading term is given by the Random Phase 
Approximation 



X3d(w,g,P±) = 



1 - 2J(u,q,p±)xid{u,q) 



Given the expression i|52|) for the dynamical susceptibil- 
ity we may obtain an estimate for the transition temper- 
ature T c below which three dimensional magnetic long 
range order develops. T c is defined as the temperature at 



which a zero frequency pole develops in X3d- Given that 
Xid(0, q) is peaked at q = and J is peaked at q = and 
p = Q we obtain the following condition fixing T c 



l-2J(0,0,Q) X id(0,0) =0. 



(53) 



Replacing J (0,0, Q) by l|47J) we arrive at the following 
equation determining the transition temperature T c 



m 



— w 2.887>d 



1 + 0A29y/Z-Z c In 



m 



(54) 



Let us consider the two limiting cases in which either 
the regular l|45(l or the singular part l|46|l of the exchange 
dominates and drives the transition. The first case occurs 
if we are very close to the point where the Fermi pockets 
are first formed and Z — Z c <C (ln(Kg))~ 2 . Then the 
transition temperature is roughly equal to 



— re 2.887 KnWln 



A 



(55) 



The second case occurs if Z — Z c ^> (ln(Ko)) 2 and then 



1.239 nl5\ 



'In 



" A " 




1 


In 


uikqS 





(56) 



where 6 = v ' Z — Z c 



V. EFFECTIVE THEORY AT LOW ENERGIES; 
THE RESIDUAL INTERACTIONS 

Now we are in position of writing the low-energy ef- 
fective action for the metallic state. This effective action 
describes interactions of the low-energy modes i.e. the 
coherent fermions and the order parameter field g a \,. In 
the previous section we calculated the interchain coupling 
for the g-field. It contains a part coming from states far 
from the chemical potential and the part with logarith- 
mic divergences coming from the states close to the Fermi 
surface. We can isolate the first piece and include it into 
the effective action of g 

s sp = J2 W ^} 

n 

+ X^ ml IdrdxTt: g m (r,x)g[(r,x) , (57) 



(52) where to first approximation 



^nl 



7 2 



7r 2 mao nl 



(58) 



This part of the action plays the role of the sigma model 
in the spin-fermion model by Pines and Chubukov (see^ 
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and references therein). Taken in isolation, this model 
has an instability at some temperature T c which can be 
estimated from the RPA expression for the dynamical 
magnetic susceptibility in complete analogy with our cal- 
culation in section llVBI Since the coherent fermions are 
low-energy excitations, they cannot be simply integrated 
out, but their interaction with g-field should be added to 
the action. The Fermi surface of the coherent fermions 
is determined by the equation 



G o (0,q) tj.(ki) = 1, 



(59) 



and consists of four pockets (two electron-like ones and 
two hole-like ones) as is shown in Fig [31 Let us consider 
the situation where the scale of the interchain hopping to 
is very slightly larger than the minimal value to required 
for the formation of a Fermi surface 



to = ta + St 



11 + 5V5 



Zo 



5t 



(60) 



The electron and hole pockets are then shallow and 
anisotropic and the Fermi surface is determined by the 
equation 



where 
E(q,k x ) 

E 

An 



Anm 



E(q,k ± )=0 



0.352<ft , 
0.543 , A ± 



vqo 
m 



1 + V5 



1.27202 



(61) 

- E , (62) 

=s 1.27202, 

(63) 



The electron pockets are formed at (kp + q, k±) and 
(—kp — q, kx) whereas the hole pockets are located at 
(kp —q, Q + kj_) and (— kp+q, Q + kjJ, where q and k^ 
are determined from 1(61(1 . Let us denote the annihilation 



operator of the coherent fermions by ^(r, q, kx)- The 
soft modes occur in the vicinity of the electron and hole 
pockets and it is convenient to decompose \I/(t, g, kx) ac- 
cordingly. We denote by R e (r, q, p±) and L e (r, q, px) the 
annihilation operator in the vicinity of the electron pock- 
ets and (q, pj_) is the deviation from (±kp, 0). Similarly 
we denote by Ru(t, q, px) and Lh(r, q, p±) the annihi- 
lation operator in the vicinity of the hole pockets and 
(q, p±) is the deviation from (±fcp, Q). 

From Ea. 1)62(1 we determine the particle density asso- 
ciated with a single pocket is 



n w 0. 027 a~ 2 q i4{Z 



\3/2 



(64) 



The liquid of quasi-particles becomes degenerate at tem- 
peratures of order of Eq. Comparing it with the transi- 
tion temperatures lfS^|) , ljSl)j l we conclude that the degen- 
erate metallic state exists only at 

Z - Z c > k 2 q . (65) 
corresponding to na\ ^> 0.027qoK^. Close to the Fermi 
surface the Green's function ((2"5)l can be approximated 
as 



G(LO,±k F + q,k ± ) 
G(uj,±k F + q,k 1 _ + Q) 



Z 2 



iuj — E(±q, kj 
Z 2 



iuj - 



where 



vqo 



0.352 . 



, (66) 
, (67) 

(68) 



The expressions l(67|l exhibit the particle-hole symmetry 
characteristic of our model at half-filling. As usual, we 
include the residue Z in the coupling constant and re- 
place the fermionic action by the action of four compo- 
nents of free fermions. The effective action describing the 
fermions is then given by 



Sf = a\J -1—3- [R* a Jr,k) (d T - £7«(k)) i? a , Q (r,k) + L* a Jr,k) (d r - E^(k)) L 0j „(r,k)] , (69) 



where k = (q, kx), a =T>Ij a = e,h and where 



„ T Ant = ik,k' Y\ R*a, a (T,k)L a ^(T,k')g a p(T,k-k') 

E*> L (k) = E{±q, kx) , E^ L (k) = -E( T q, k ± ). (70) Q t^ 

+ Iw,Q+k'R*e, a {T, k)L M (r, k')g af3 {T, k - Q - k') 

The fermion-spin vertex is described by the action + / Q+k k , R* ha (r, k)L e ^(r, k')g af3 (r, k + Q - k') 

+ h.c. , 

, f dr d 3 k d 3 k' „ ,„ s v z Zn 

Sint = a± J (2^)« int ' ( } /k < k ' = 2 ^^^( k )^( k ') • ( 72 ) 
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All wave vectors in the above formulae lie close to the 
non-interacting Fermi surface and therefore their longi- 
tudinal components are small in comparison to it: \q\ <C 
■na,{). The entire approach is valid only when the volume 
inside of the Fermi surfaces is small. One then can neglect 
the momentum dependence of the exchange constant in 
Ea.(|72p. The sign of the exchange constant depends on 



the "pocket index" a, b 



7m 



-1 1 



(73) 



where 7 is a constant. The interaction can be cast in the 
form 



J 



7m 



K, Q (k)L 0)/3 (k')<Mk ~ k') - {i?; a (k)L fel/3 (k') + ^ a (k)L eJ3 (k')}g aP (k - k' - Q)j + h.c. (74) 



The value of the coupling constant 7 can be extracted 
from Eqns (|42[1 and (|7(j|l by noting that it is this inter- 
action which gives rise to the logarithmic singularity in 

J(Q) 



J sing (0,0,Q) 



9 2 2 P(0) 

-27 TO 

a 



hi 



(75) 



Here p(0) is the density of states per species at the Fermi 
surface of coherent fermions 



p(0) = lim ao 

LU— >0 



0.539 



d 3 k 

(2tt) 3 



a k 



-ImG(w, ftp + q, kjj 

7T 



(76) 



St = \ drd 3 x *t ( T) x ) <J (/ ® JJSr + {I®t z ) 

Shit 



The result is that close to the transition we have 



TO 



(77) 



where t 3> to is the hopping along the chains and the 
constant of proportionality is of order 1. Though 7 is 
never small, the small parameter Kq appears every time 
one integrates over the transverse momentum. Hence 
the magnitude of 7 is not a problem. The effective action 
describing the metallic side of the Mott-insulator to metal 
transition is given by Eqns 1)57(1 . (|69|l and (|74|l . We find 
it instructive to write it down also in position space 



En 



if. 



vi 



2Mn 2Mi 



^ /" drd 3 x *t (t,x) ({ T + ^ [exp(-2z go xr 2 ) - r 1 expHQ • x ± )]} <? Q0 (r,x) + h.c.) 4>(t,x) 



(78) 



Here we have taken the continuum limit in the directions 
perpendicular to the chains and introduced a field $J = 

{^* R ta^Rt,oi^ L t^4>* L t ta ), where = exp(«g a;)- 
We employ a tensor-product notation, where the first 
space is associated with the "right/left" index and the 
second space with the "e/h" index. The Fermi surfaces 
of electrons and holes are shifted to the origin and su- 
perimposed. The spin action S[g] is given by Ea. (|57ll . 
Alternatively, one may use the Abelian representation 
given by Eq.@, with g defined by (fT3f) . 



A. Marginal Fermi liquid 

As we shall now demonstrate, at temperatures higher 
than the Neel temperature T c this metal is, in fact 



a Marginal Fermi Liquid^. The following discussion 
closely parallels the analysis given by Chubukov et al for 
the spin-fermion model (see, for example Rcf. 2]). Let 
us consider diagrams for the Green's function of right 
moving electrons. We expand around uncoupled chains 
and take both the spin-fermion coupling and interchain 
spin-spin exchange into account perturbatively. The ele- 
ments of the diagram technique for the fermionic degrees 
of freedom are as usual, whereas the building blocks in 
the spin sector are the connected 2n-point spin correla- 
tors for a single chain. In diagrams that do not contain 
closed fermionic loops or the interchain exchange such 
as the ones in Figs 7a and 7b, the spin correlations are 
independent of the transverse wave vector. This means 
that each fermion Green's function is integrated over kj_ . 
This integral does not differ significantly from the inte- 
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gral over all momenta and as a result is independent of 
qu , corresponding to a Green's function that is local in 
real space: 



dk 2 , 



(27T) 



const i«Q sgn(w) . (79) 



As (|79|) is independent of q, we may integrate the spin 
correlator in the diagram of Fig. 7a over q. This makes 
the spin correlator local. As a result the contributions 
to the self energy which do not contain closed fermionic 
loops or interchain spin exchange are approximately mo- 
mentum independent. Then the self energy calculation 
becomes essentially a local problem like the problem of 
electron-phohon interactions in metals and superconduc- 
tors (the Eliashberg theory In fact, such an approach 
works under less stringent conditions, namely, when the 
spin excitations in the transverse direction are much 
slower than the quasi-particles. Therefore the diagrams 
generating a kj_ dependence of the spin-spin correlators, 
such as the ones in Fig. 7c) and Fig. 7d) do not af- 
fect the result for the electron self energy even close to 
the transition. Once such diagrams are neglected, we get 
an expansion where a factor Kq is associated with each 
fermionic line (originating from the integration over kj_, 
as in Ea. H79(l . Since £ depends only on frequency, mak- 
ing these lines fat does not change the result (|79|l and no 
self-consistency is required. The contribution from the 



The contribution of the diagram in Fig. 7a) to the self- 
energy is then 



max{tj, T} 



LO) CX 



iKnW In 



dr 



T 



t\ sin(7rTr)| 



In 



(81) 



Here E 1 serves as a short-time cutoff in all integrals. 
The diagram in Fig 7b) involves cumulants of the type 



(82) 



and gives a contribution 



OC Kq 



T13T24 



T12T14T23734 



T12T23T34 

1 



T12T34 



F14T23I 



oc u? (83) 



to the self energy. Various other contributions are zero 
because some local cumulants vanish 

(( e f*(Ti) e -ie(r 2 ) e |e(r 3 ) e -i*(T4)^ = ( 84 ) 

Equation 4point shows that contributions from higher 
cumulants can be neglected at small frequencies. As a 
result, the only essential contribution to the self energy 
comes from the diagram Fig. 7a and is given by Eq. 18111 . 



VI. 



THE ORDERED STATE 



a) — 

t 



"O- 

\ 



b) — i° v — 



9 9 9 9 

c) > — d) > — 



FIG. 7: Diagrams for the quasi-particle self energy of right- 
moving electrons. The lines with arrows represent the 
fermionic Green's functions of right and left moving electrons 
and holes. The 2n-point vertices denote cumulants of the 
matrix fields g and g^ . 



diagram Fig. 7a contains the correlation function 



©On) e -|e(T 2 ) 



(in [A]) 



1/2 



j sin(7rT[ri - r 2 ] 



(80) 



As we have discussed in the previous sections, the 
system undergoes an antiferromagnetic transition at a 
temperature much smaller than the Mott-Hubbard gap 
T c ~ ttikq. Once T c becomes small compared to the 
quasi-particle Fermi energy Eq, one can distinguish be- 
tween metallic and insulating behavior. As we have 
demonstrated, the corresponding metal is rather unusual, 
being in fact a Marginal Fermi Liquid. Below T c however, 
the system becomes either an insulator (for zero doping) 
or an ordinary Fermi liquid. Indeed, at zero doping the 
electron and hole Fermi surfaces are nested and the or- 
dering occurs at the antiferromagnetic wave vector in the 
transverse direction (recall that the chains run along the 
z-axis) such that 



(Si.a/sO)) = a af} ■ M(-l)'- +, » 



(85) 



Here the components of a are the Pauli matrices and M 
is the ordering vector. In the mean-field approximation 
the fermionic spectrum is gapped 



(^(k)) 2 + 7 2 m 2 |Mf 



(86) 



At non-zero doping our approach still holds provided 
the chemical potential lies inside the Mott-Hubbard gap. 
There is no nesting any longer and the magnetic ordering 
does not open a gap in the quasi-particle spectrum. As 
usual, magnetic fluctuations interact with quasi-particles 
through gradient vertices and these interactions are weak. 
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VII. CONCLUSIONS 

The main result of this paper is a formulation of a 
self-consistent description of the hybrid state of 3D quasi- 
particles interacting with magnetic collective modes. The 
derivation is done for a toy model of half filled Hub- 
bard chains weakly coupled through a long range inter- 
chain hopping. A certain artificiality of the model was 
necessary to ensure the self-consistency of our approach 
through the presence of a small parameter Kq. We have 
also neglected the long-range component of the Coulomb 
interaction, which plays an important role in determining 
the character of the metal-insulator transition. In reality 
a long-range interaction may lead to an instability of the 
small FS phase, though for small Mott-Hubbard gap its 
influence is diminished by the presence of a large dielec- 
tric constant. In this case the first order MI transition 
line may terminate below the antiferromagnetic transi- 
tion line of Fig. 1. 

The resulting low energy effective theory is of the 
Eliashberg type: the interaction between quasi-particles 
and collective modes leads to strong retardation effects 
resulting in a strongly frequency dependent quasi-particle 
self energy. In the present model this takes place not 
just at the 'hot spots', as in the spin-fcrmion model of 
Chubukov and Pines, but on the entire quasi-particle FS. 
This makes the present model a candidate for the de- 
scription of 'bad' metals. The fact that in our model the 
electron self energy is of the Marginal Fermi Liquid form 



is not universal and is determined by the particular spin 
fluctuation spectrum. 

As we discussed in the introduction, our theory pro- 
vides an example of a state where the number of carriers 
is unrelated to the volume of the FS. Though this idea 
is well established (see e.g. the textbook |5J), its micro- 
scopic realization was restricted to superconductors (the 
example given in Ref.j^l). Our model provides another 
example. It also demonstrates that one does not need ex- 
otic ground states to have a small FS, as was suggested 
in Refs [23], [HI- The small FS phenomenology can be 
generalized beyond our model. In general there is no 
a priori reason for the Fermi surface even to be closed; 
for instance, Ref.|2^] describes a state with a truncated 
Fermi surface observed in ARPES experiments on un- 
doped cuprates^. 
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